Thermal spin fluctuations in spinor Bose-Einstein condensates 
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We study the thermal activation of spin fluctuations in dynamically-stable spinor Bose-Einstein 
condensates. We analyze the specific cases of a non-dipolar spin-1 condensate in m — 0, where 
thermal activation results from spin-changing collisions, and of a Chromium condensate in the 
maximally stretched state m = —3, where thermal spin fluctuations are due to dipole-induced spin- 
relaxation. In both cases, we show that the low energy associated to the spinor physics may be 
employed for thermometry purposes down to extremely low temperatures, typically impossible to 
measure in BECs with usual thermometric techniques. Moreover, the peculiar dependence of the 
system's entropy with the applied Zeeman energy opens a possible route for adiabatic cooling. 
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I. INTRODUCTION 

Ultra-cold spinor gases, formed by atoms with multiple 
internal states, have attracted a major interest due to the 
rich physics resulting from the interplay between internal 
and external degrees of freedom [l[. In particular, spinor 
Bose-Einstein condensates (BECs) have been throughly 
investigated forgases of different spins, including spin- 
1 (F = 1 87 Rb and 23 Na @), spin-2 (F = 2 87 Rb Q), 
and spin-3 ( 52 Cr @,@). 

Especially interesting is the possibility of a coherent 
population transfer between the different internal compo- 
nents. This is accomplished through the so-called spin- 
changing collisions, which, although conserving the to- 
tal spin projection of the colliding pair of atoms, in- 
duce a redistribution among the various internal lev- 
els These collisions originate from the differences 
between intra- and inter-component scattering lengths. 
Since these differences are very small in usual experi- 
ments, spin-changing collisions are hence associated to 
very low energy scales, typically much smaller than the 
overall chemical potential. 

The case of Chromium deserves a separate discus- 
sion. Chromium is not only an example of a spin-3 gas, 
but presents as well large magnetic dipole-dipole inter- 
actions (DDI) 0- Contrary to the usual contact-like 
isotropic interactions which as mentioned above preserve 
spin projections, DDI arc anisotropic and hence allow for 
spin relaxation, a phenomenon which has been explored 
in Chromium BECs in recent years 0, Q . As for the case 
of spin-changing collisions in non-dipolar spinor BECs, 
spin relaxation in Chromium is associated as well to very 
low energy scales. 

In this paper, we study spin fluctuations in spinor 
BECs, and in particular the dependence of these fluc- 
tuations on temperature. BEC thermometry is typically 
performed using timc-of-flight measurements, either from 
the expansion velocity of the thermal cloud or from a 
bimodal fitting which allows to establish the ratio be- 
tween condensate and thermal cloud. These techniques 



fail however for low temperatures for which the thermal 
population is small compared to the number of particles 
in the condensate. The use of low-energy phase fluctu- 
ations for thermometry in a two-well scalar condensate 
was pioneered in Ref. In this paper, we show that the 
very low energy associated to both spin-changing colli- 
sions and spin-relaxation opens as well interesting possi- 
bilities for thermometry purposes down to extremely low 
temperatures. Moreover, we show that the dependence 
of the entropy of the gas on the Zeeman energy allows 
for a possible mechanism for adiabatic cooling. 

Our paper is structured as follows. In Sec. [TT] we an- 
alyze the case of a stable spin-1 BEC prepared in the 
m = Zeeman sublevel. We study by means of the cor- 
responding Bogoliubov analysis the thermally activated 
spin fluctuations resulting from spin-changing collisions. 
Sec. IIIII is devoted to the case of a Chromium conden- 
sate prepared in the maximally stretched Zeeman state, 
m = —3. This case differs significantly from the spin- 
1 case, since the thermally activated spin fluctuations 
result from the spin-relaxation induced by the DDI. We 
show that in both scenarios spin fluctuations may be em- 
ployed for deep temperature thermometry and adiabatic 
cooling. In Sec. IIVI we summarize our conclusions. 



II. SPIN-1 CONDENSATE 

A spin-1 BEC is described by the Hamiltonian 0] 
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where ^ m (r)(^] n (r)) is the annihilation (creation) oper- 
ator of boson in the Zeeman state m = 0, ±1 at position 
r, M is the atomic mass, Vt rap (r) is the external trap- 
ping potential, and F = (F Xl F y , F z ) is the spin vector 
operator (F x ,y,z are the spin-1 matrices). The couplings 
c = 4nh 2 (a '+ 2a 2 )/3M and c 2 = 4-Kh 2 (a 2 - a )/3M, 
are expressed in terms of the scattering lcgths 0,0.2 which 
characterize low energy collisions with total spin and 2, 
respectively. The term qm 2 denotes the quadratic Zee- 
man shift in an external magnetic field B , with q = 
/i^i?o/8Chf s , with Chfs the hypcrfinc coupling strength. 
Note that collisions conserve the total spin projection and 
hence the linear Zeeman energy is a conserved quantity 
which may be gauged out. 

We assume a stable m = BEC (we will discuss 
the conditions for stability below). Assuming small 
fluctuations around the condensate solution, the spinor 

BEC is described by the field operator ^ = (0,tp o ,0) + 
(5^-i, 5^0, dty+i), where the BEC wavefunction ipo ful- 
fills the Gross-Pitaevskii equation 
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^o(r) = M)(r) , (2) 



where fi is the chemical potential, and no(r) = \ipo(r)\ 2 . 
Retaining up to second order in the fluctuations, we ob- 
tain an effective Hamiltonian for S^>±i: 

Hi = ! dr S^ m H cS 5^ m 

m=±l 

+ c 2 dr n + h.c. 



(3) 

where H cS = -h 2 V 2 /{2M) + V(r) + q, with V(r) = 
Vtrap(r) + (c + C2)n (r) — (x. In homogeneous space, 
Vtrap = 0, no is constant and we may move to momen- 
tum (k) space where the Hamiltonian becomes: 

H x = f dke k ]T 5& m (k)6if m (k) 

•> m=±l 

,J dk [<5* +1 (k)5*_ 1 (-k)+h.c. 
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C2 n o + 1- Using the symmet- 
ric and antisymmetric bosonic operators Su_, = 
f t^+i (k) ± 5^_i(k)V and applying a Bogoliubov 



transformation Sk,^4k = f^B^ 
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where E k = \J e k — c 2 n o * s * ne spectrum of s P m excita- 
tions. Note that the m = BEC is stable as long as E k 
is real, which demands q > q cr , with q cr = (\c 2 \ — c 2 )riQ. 
In equilibrium at a temperature T (/3 = 1 /ksT with fee 

the Boltzmann constant) (B^B^) = (e~ pE ^ - 1 
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FIG. 1. Ratio N±i/Nq in the trap as a function of q/q cr and 
ksT/fi for the case discussed in the text. 



with i = S,A, and (B^B^ i] ) = 0. Re-expressing 
5i5f±i as a function of the quasi-particle excitations, we 
obtain the density in m = ±1 

1 



n±i(T) 



(2tt)= 



dk n±i(k, T) . 



with n±i(k,T) = (5^] tl 6^± 1 ) given by 



n±i(k,T) 
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the same for both components. 

We consider at this point a trapped BEC in the 
Thomas- Fermi (TF) regime with density profile no(r) = 

£ (1 - rl/Rl - z 2 /R 2 z ), where R ±iZ = ^/(M^~) 

are the TF radia. For a sufficiently smooth density 
profile, we may employ the local density approxima- 
tion (LDA) , associating to each value of the local density 
no(r) the corresponding excitation spectrum for the ho- 
mogeneous case with that density, E k (r). We may then 
evaluate the corresponding local density n±i(r) from the 
expressions obtained above. The total number of atoms 
inm = ±l is obtained by integrating their local occupa- 
tion over the density profile of the trap: 



(JV±i)= / d 3 r(n±i(r)>, 



(8) 



The critical value of the magnetic field, q cr is calculated 
at the trap center. 

We consider in the following the specific case of F = 1 
87 Rb, for which ao = 101. 8as and a 2 = 100. Aclb (with as 
the Bohr radius). As a consequence, c 2 = —4.6 x 10 _3 Co 
is small and negative, which provides a critical q cr //J, = 
2 1 C2 1 /co = 9.25 x 10~ 3 . For simplicity we consider that 
the atoms are confined in a spherically symmetric trap, 
with a harmonic frequency u> = 2tt x 50Hz. For a typical 
value of N = 10 5 atoms, the density at the trap center 
becomes 10 14 cm -3 . 

In Fig. [TJ we plot the total number of atoms (N±i)/Nq 
in the trap as a function of ksT/n for different values 
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FIG. 2. Entropy as a function of ksT f [i and q/q cr for the 
same case of Fig. [T] The black curves indicate various iso- 
entropic curves, the value of S/ks is indicated. 



j-2i2/fc2i2 \ 

I + 2cono J , corresponding to S^q fluctuations, 

also contribute in principle to the system entropy. How- 
ever, if C2 -C Co (as it is typically the case), the entropy 
contribution of the scalar excitations may be neglected 
for ksT <C fJ-- In particular, we have checked that for q in 
the vicinity of q c the entropy contribution of the scalar 
modes is negligible for ksT/fi < 0.25. For a trapped 
gas, we employ again local density arguments, calculat- 
ing the local entropy density associated to each position 
s(r). Fig. [2] shows curves of equal entropy S — J dr s(r) 
as a function of temperature and q/q cr for the same case 
discussed in Fig. [TJ Interestingly, the isotropic curves 
bend to lower temperatures when approaching q cr . As 
a result, an adiabatic variation of the applied magnetic 
field may allow for adiabatic cooling. For example, in 
the figure, starting with q = 5q cr at T = OAfx/ks, the 
system may decrease its temperature down to 0.05/i/fcs 
when approaching q cr . 



III. CHROMIUM CONDENSATE 



of q/q cr - As expected, when approaching the critical q cr 
the spin population is enhanced at low temperatures, due 
to the very low energy associated to spin excitations. At 
q = q cr the population of m = ±1 (always for N = 10 5 
atoms in the BEC) is larger than 10 atoms for tem- 
peratures larger than 0.02/i/kB- Although these num- 
bers may be very small to be observed experimentally, a 
simple procedure may significantly enhance the experi- 
mental resolution of such small populations. If the sys- 
tem is abruptly brought into instability by sweeping into 
q < q cr , spin excitations grow in a process similar to para- 
metric amplification in non- linear optics (lol | . Within the 
so-called linear regime, spin fluctuations are hence expo- 
nentially amplified with a growth rate T. As a result 
the population in ±1, is enhanced in time in the form 
N±i(t) ~ iV±i(0)e rt , where iV ±1 (0) is the population in 
m = ±1 prior to the destabilization, i.e. that depicted in 
Fig.[TJ A subsequent Stern-Gerlach arrangement in time- 
of- flight , allows for a separate imaging of the different 
Zeeman components. As a result even very small N±i(0) 
may be experimentally resolved, opening the possibility 
of employing thermally activated spin fluctuations as a 
thermometer down to temperatures close to 0.01/i/fcs 
where other thermometry methods typically fail. 

Interestingly, the dependence of spin excitations on q 
may be employed as a possible mechanism for adiabatic 
cooling. For an homogeneous gas, the entropy density 
of the system is readily calculated from the spectrum of 
elementary excitations 
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We consider in the following the case of a Chromium 
BEC, which constitutes an example of spin-3 gas in 
which, crucially, strong magnetic dipolc-dipolc interac- 
tions (DDI) induce spin relaxation. As a consequence the 
low temperature physics is considerably different com- 
pared to that of spin-1 BECs discussed in the previous 
section. For the spin-3 Chromium BEC the Hamiltonian 
acquires the form H = Hq + V sr + Vdd, where 



dr^ m 
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where the linear Zeeman energy is included since, con- 
trary to the spin-1 case, the global spin projection is 
not conserved. Note also that in principle Chromium 
presents no quadratic Zeeman effect due to the absence of 
hyperfine structure. However, an effective quadratic Zee- 
man effect may be induced by optical or microwave dress- 
ing [ll|, [l2j ■ Contact interactions are described by [l3| : 



V. T = 



dr 



c n 2 (r) + Cl F 2 (r) 
+c 2 P Q (r)+c 3 6 2 (v) 



(11) 



Note that scalar excitations of the m = condensate, 



where the symbol :: denotes normal order; n(r) = 

EmAW^W is the total density, F\r) = 

Em,n&W^»W. F 2 (r) = £,:(^(r)) 2 ' V a (r) = 

7E ra J"l) m+ "4(r)^,(r)iWin(r), 6«(r) = 

E m , n ^)(S^S=) m J n (r), 6 2 (r) = V ( ! (a. ; ;r:) 2 . 
The interaction constants are cq = (— II92 + 8I54 + 
7ffe)/77, ci = (g e - 92) /18, c 2 = go + (-55g 2 + 27g A - 
5<7 6 )/33 and c 3 = g 2 /126—g i /77+g 6 /198, where as in the 
previous section, g$ = Anh 2 as /M for the channel of total 
spin S. For 52 Cr, a 6 = 112ob, c = 0.65g 6 , c\ = 0.059g 6 , 
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c 2 = .90 + 0.374(76 and C3 = — 0.002g 6 (the value of ao is as 
yet unknown). Finally, the DDI acquires the form (l3l |: 



To Cd 



dr dr' 

r _ r'|3 



: J- Z ,,(r,r')r 20 (r - r') + ^.-(r, r')F 2 i(r - r') 
+J' z , + (r,r')r 2 _ 1 (^P)+ J-_,_(r,r')r 22 (r^?) 



+J-+,+(r,r')r 2 - 2 (r-r') 



where Y2 m {r — r') are the spherical harmonics, 



2 r 



(12) 



•F„,,(r,r') = y-[3F z (r)^(r') - F(r) • F(r') 

J- 2 ,±(r,r') = ±[F±(r)F,(r') + A(r')F±(r) 
^±,±(r,r') = F±(r) J F±(r'), 



(13) 



and j±(r) = F x (r) ± iF y (r). The DDI arc characterized 
by the coupling constant = \XQir B g\l Al-k , where jtto is 
the vacuum magnetic permeability, /^b the Bohr mag- 
neton and gL the Lande factor. For 52 Cr, g L = 2 and 
c ( ; = 0.00456- Note that the DDI do not conserve spin 
and orbital angular momentum separately allowing for 
spin relaxation processes, in which e.g. atoms in m = —3 
arc transferred into m = —2 

In the following we consider that the linear and 
quadratic Zeeman effects arc chosen in such a way that 
only the two lowest states of the Zeeman manifold, m = 
—3 and m — — 2 contribute, whereas spin relaxation to 
other m states is energetically suppressed. In this simpli- 
fied scenario, we assume a condensate of m = — 3 atoms 
with small spin fluctuations populating the m = — 2 
component. This system can be described by the field 
tp ~ ?/)_3 + <5?/>_3 + 8tj)—2, where the BEC wavefunction 
fulfills the GP equation 
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gn-3(r)ip-a(T) 



+ Vtrap(r) -3p+9q 



^-s(r) 
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with g — cq + 9ci + 8IC3. 

In homogeneous space, Vt ra p( r ) = 0, ips(r) = 
and fi = 9q — 3p + gri-3. Moving into momentum space, 
5ip m (k), we introduce the operators 

d_ 3 ,±(k) = ^_3(k)±^I 3 (-k) 

0- 2 ,±(k) = c-^ fc ^_ 2 (k) ± ^IjC-kJe*** , (15) 

where we have introduced spherical coordinates k = 
(k, 9 k , 4>k)- These operators are governed by a set of cou- 



pled Hciscnbcrg equations: 
ih 




where 




0-s,-(k) 
0-2,-(k) 

0-3,+ (k) 
0-2,+(k) 



M 2 i M 2 2 



+ 2gn-z + 247rc (i n_ 3 (3 cos^ 9 k - 1) 
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367r\/-c d n_ 3 sm6» fc cos9 k , 



21,2 



M 2 



2M 



j-V- 47rc d n_ 3 (l - 3 sin 2 9 k ) 



(19) 



Note that contrary to the spin-1 case, scalar fluctuations 
(given by 8^-3) couple with spin fluctuations (given by 
8^-2) at first order. As a consequence, the elementary 
excitations have a hybrid scalar/spin character absent in 
the spin-1 case. 

The corresponding Bogoliubov excitations may be 
written as a linear combination of the operators above: 



/ A+(k) \ 
A_(k) 
A + (-k) 

VAL(-k); 



= T(k) 



/6-s,-(k)\ 
0-2,-(k) 
0- 3 ,+ (k) 



(20) 



where 



and 



T(k) = 

a 2 ± (k) 

/3±(k) 
7±(k) 
5±(k) 



V6_ 2 ,+(k) / 

/a+(k) /?+(k) 7+ (k) * + (k) 

a_(k) /3_(k) 7 _(k) 5_(k) 

a+(k) /?+(k) - 7+ (k) -<J+(k) 

\a_(k) /3_(k) - 7 _(k) -$_(k) 



(21) 



B 2 (k)R n (k)R 2 2(k) 



a±(k) 
a±(k) 
/?±(k) 



4S±(k) 
££(k)-A(k)) 2 i?n(fc) 
£ 2 (k)-A(k) 



B 2 (k)i? 22 (fc) 
1 



B(k) 

£±(k) 
i?n(fc) 

g±(k) 
i? 22 (fc) ' 



(22) 
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The excitation Hamiltonian acquires then the form 

H = f dk J2 E K (k)ki (k)A K (k) . (23) 



with eigenenergies: 

E 2 ± (k) = ±(A(k)+D(k)) 



±^(A(k)-D(k)) 2 +4C(k)B(k), (24) 

where A(k) = R n (k)M n Qz), B(k) = R ll (k)M 2 i{k), 
C(k) = i? 22 (fc)Af 12 (k), and D{k) = R 22 {k)M 22 (k). The 
modes are stable, i.e. possess real eigen-energies, if 
U < U cr = —Prized/ g§- Note that stability is just gov- 
erned by the interplay between Zeeman energies and the 
spin relaxation due to the DDI. 
Defining the matrices: 



A 



r 



1 



1 







— a_ 


a + 


8- 




-7- 


7+ / 



7+<5- — 7_<5_ 

we may express the population in m = —2 in the form 
A 2 + M- 



(25) 



(n_ 2 (Mk)> = 

gfc B T _ \ & k 

where we introduce the amplitudes 



- 1 



^2+(k) = - 

4>-(k) = I 



r^kJ + A^k) 
rL(k) + A^ 2 (k) 



Q-2 (26) 

(27) 
(28) 



and the zero temperature quantum fluctuations 

n 2 



Q- 2 (k) 



r 21 (k) + A 21 (k) 



r 22 (k) + A 22 (k) 



(29) 



Contrary to the spin-1 case discussed in the previ- 
ous section, (n_ 2 ) has in general an angular depen- 
dence, which results from the anisotropy of the DDI. 
We may quantify the anisotropy of the spin population 
by means of x = J(3cos 2 #fc — l)(h- 2 (k,6k))d 3 k. An 
isotropic distribution is characterized by \ — 0, whereas 
positive values indicate a distribution preferentially ori- 
ented along 6k = 7r/2. The anisotropy x presents an 
interesting dependence as a function of temperature and 
U/U cr , depicted in Fig. [3] for the case of n_ 3 = 10 14 
cm~ 3 , U cr /fio = —0.05. At low T, x nas small pos- 
itive values, increasing when U/U cr increases (for very 
low T < O-Ol/i/fce, x acquires a maximum for interme- 
diate U values). For larger T, x < for low U/U cr indi- 
cating a momentum distribution oriented along 6k = 0, 
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FIG. 3. Anisotropic x as a function of fesT/^i, 
U/Ucr, for an homogeneous 52 Cr BEC with n_3 = 



and 
10 14 



particles/cm . The black lines describe configurations with 
the same anisotropy, the corresponding x values are indicated. 



whereas at larger U /U cr the distribution becomes basi- 
cally isotropic. 

We consider at this point the case of a trapped 
Chromium condensate. In general, LDA must be care- 
fully considered, due to the long range character of the 
DDI. However, when the characteristic length of this 
interaction, add = 3Mcd/fi, is much smaller than the 
typical length of the condensate harmonic trap, a^. = 
y/h/Muj, the LDA can still be used to calculate the num- 
ber of particles in the m = —2 state, as long as the den- 
sity profile of the m = —3 BEC varies smoothly with r. 
This approximation allows to estimate the total number 
of atoms in m = — 2, but is of course not appropriate to 
study its angular distribution. 

For an axially symmetric harmonic potential 
Vtrap(r) = Mu]jr 2 x + A 2 z 2 )/2, where A = uj z /uj±_ 
is the trap anisotropy. The density profile is, 

IhN 



n(r±, 



1 



R\ 



z 

B 2 



" 8ttR 2 ± R 

where k = R± / R z is found by solving [14j, [1 
c 2 V 2 



(30) 
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A 2 " 



3e 
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2e 



dd 



1 



E d d - 1 ,(31) 



with edd = 127r 2 Cd/g, and 

1 + 2k 2 3k 2 



/(«) = 



1 



(1 - K 2 ) 3 / 2 



tanh -1 yjl - k 2 . (32) 



Normalizing the total density to the total number of 
particles, one finds: 



R± = 



1/5 



(33) 
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FIG. 4. Number of Chromium atoms in m = —2 as a function 
of ksT I \x for different values of U/U cr , for the case discussed 
in the text. 



where a = \hgN / A-kMuj]_. The chemical potential is /j, = 

gn(0,0)^l — At each small volume, the local 

chemical potential is /x(r) = ji — Vtrap(r), which may be 
well approximated for Chromium by fj,(r, z) ~ gri-^ir, z), 
as for the homogeneous BEC. 

Fig. |4] shows A^_2 as a function of temperature for 
several values of the U/U cr , for the specific case of a 
52 Cr BEC with N = 10 5 in m = — 3 in a spherical trap 
of frequency u = 2n x 50Hz, where U cr is determined 
by the central density n(0,0) = 5.33 • 10 13 cm~ 3 . Note 
that add = 0.894nm <C clh.o. = 1.97/im, and hence well 
within the limits of the LDA. Close to U cr populations 
of N-2 = 10 may be attained below 0.1/i/ ks- Hence, 
as for the spin-1 case, Fig. |4] shows clearly that one may 
employ the population in m = —2 (combined with an 
abrupt jump into instability, as discussed for spin-1) for 



thermometry purposes. Finally, we should note that as 
for the case of spin-1 isotropic curves bend towards lower 
T when approaching U cr , and hence also for Chromium 
an adiabatic reduction of U/U cr may allow for an inter- 
esting cooling mechanism. 



IV. CONCLUSIONS 

We have analyzed the thermal activation of spin ex- 
citations in spinor condensates. For the case of spin- 
1 condensates, an stable m — condensate presents a 
non- negligible thermally activated population of m = ±1 
due to spin-changing collisions. For the case of a sta- 
ble Chromium BEC in m = —3 dipole-induced spin- 
relaxation leads as well to thermal activation, which 
contrary to the spin-1 case, acquires an intriguing 
temperature-dependent anisotropy. For both cases we 
have shown that the spin population may be employed 
at very low T <C /i/fcs as a possible mechanism for deep- 
temperature thermometry. We have shown as well that 
an external adiabatic variation of the Zeeman energy may 
be employed to achieve an adiabatic cooling mechanism. 
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